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Abstract 

We explore the temperature dependence of the heavy-quarkonium interaction based 
on the Bhanot - Peskin leading order perturbative QCD analysis. The Wilson co- 
efficients are computed solving the Schrodinger equation in a screened Coulomb 
heavy-quark potential. The inverse Mellin transform of the Wilson coefficients then 
allows for the computation of the IS and 2S heavy-quarkonium gluon and pion to- 
tal cross section at finite screening/temperature. As a phenomenological illustration, 
the temperature dependence of the IS charmonium thermal width is determined and 
compared to recent lattice QCD results. 
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1 Introduction 



The Debye screening between two opposite color charges is clearly seen in the 
QCD static potential computed at finite temperature T on the lattice [jj. Con- 
sequently, heavy quark bound states (which we call <3?) may no longer exist well 
above the deconfinement critical temperature T c , of order 200—300 MeV 0. 
This has made the heavy-quarkonium suppression in high energy heavy-ion 
collisions (as compared to proton-proton scattering) one of the most popu- 
lar signatures for quark-gluon plasma formation |3j,|4j. On the experimental 
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side, a lot of excitement came out a few years ago after the NA50 collabora- 
tion reported a so-called "anomalous" suppression in the J/ip channel in the 
most central lead-lead collisions (y/s ~ 17 GeV) at the CERN SPS fl- At 
RHIC energy (y/s = 200 GeV), J/ip production has been measured recently 
by the PHENIX collaboration although the presently too large statistical and 
systematic error _bars prevent one from concluding anything yet quantitative 
from these data 



The NA50 measurements triggered an intense theoretical activity and subse- 
quently a longstanding debate on the origin of the observed J / ip suppression. 
However, it became unfortunately rapidly clear that no definite conclusion 
could be drawn as long as theoretical uncertainties exceed by far that of the 
high statistics data. Indeed, both the realistic description of the space-time 
evolution of the hot and dense medium as well as the interaction of heavy- 
quarkonia with the relevant degrees of freedom (let them be pions or gluons) 
are required to be known. While the former can be constrained by global ob- 
servables, the latter needs to be computed theoretically. Several approaches 
have been suggested to determine heavy-quarkonium total cross sections, from 
meson exchange 0] or constituent quark models to the perturbative frame- 
work developed by Bhanot and Peskin 0, 0] upon which the present paper 
relies. Let us remark in particular that many recent phenomenological appli- 
cations have used the latter perturbative $ - gluon cross section to estimate 



the heavy-quarkonium dissociation or formation in heavy- ion collisions 11 



However, although derived from first principles in QCD perturbation theory, 
the Bhanot - Peskin result describes the interaction of Coulombic bound states, 
that is for which the heavy-quark potential is well approximated by the per- 
turbative one- gluon exchange potential. As indicated from spectroscopic stud- 
ies |l2| , this may be too crude an assumption to describe bound states in the 
charm or (even) the bottom sector. Furthermore, it does not take into account 
the possible effects of the medium on the heavy-quarkonium interaction. It is 
the aim of this Letter to explore how the $ interaction with gluons and pions 
gets modified at finite temperature. The paper is organized as follows. The 
general framework is first briefly recalled in Section 2. Our results are then 
detailed in Section 3 while Section 4 is devoted to a concluding discussion. 



2 Heavy-quarkonium interaction in QCD 



2. 1 Resummation of the leading-twist forward scattering amplitude 



At leading- twist, the forward heavy-quarkonium ($) - hadron (h) scatter- 
ing amplitude M<s>h is an operator product expansion of perturbative Wilson 
coefficients d 2 k evaluated in the heavy-quarkonium state and computable in 



2 



perturbation theory times non-perturbative matrix elements in the hadron 
state. It reads 



fg 2 N c \ 



where ao and e stand, respectively, for the Bohr radius and the binding energy 
for the $ system, g the QCD coupling and N c the number of colors. Each of 
the matrix elements (h\ . . . \h) in Eq. (1) is proportional to a traceless fully 
symmetric rank 2k tensor in the spin-averaged hadron state 

n w "' M2fc (p) = jf 1 ... p^ k - trace terms 

where p M is the hadron momentum. The trace terms correspond to target mass 
corrections 0{m 2 l /e 2 ) which are neglected here as we shall deal only with 
pions in the present approach. Note that such corrections were systematically 
included in Refs. Il4j and proved relevant only slightly above the threshold 



for the quarkonium-hadron interaction process. The matrix elements can be 
written as 

(h\ l -F^{iD G ) 2k ~ 2 F v Q \h) = A 2k Il - (p) = A 2k \ 2k (2) 

where A = p° is the hadron energy in the $ rest frame and the A 2 k coefficients 
are the Mellin transform of the unpolarized gluon density G h in the hadron 
target 0,0 

A 2k = f'-x^G^x). 
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Plugging (2) in (1), the leading-twist forward scattering amplitude can be 
written as 

M* h (\) = (^) a 2 Q e E d 2k A 2k (X/e) 2k . (3) 
\ 167r / k >i 

Expressing the Wilson coefficients in terms of their Mellin moments, 

- 1 dx 



d 2k = [ — x 2k d{x), 
Jo X 



the power series (3) can be conveniently resumed and continued analytically 



throughout the whole complex plane of energies |1J]. This allows for the com- 
putation of the imaginary part of the forward scattering amplitude along the 
real axis, A > e, 

^M W =(m^j^^)s(±). ( 4) 
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Dividing Eq. (4) by the flux factor A leads to the total heavy- quarkonium cross 
section via the optical theorem 

1 r 1 

c*fi,(A) = — Im A4.[X) = / dx G{x) a^ g (x\), (5) 
A jo 

where the heavy-quarkonium gluon cross section is defined as 
with the gluon energy uj = A x in the $ rest frame. 

The Wilson coefficients need first to be computed in an arbitrary heavy-quark 
potential and later be inverse Mellin transformed in order to determine the 
heavy-quarkonium gluon Eq. (6) and hence the heavy-quarkonium hadron 
Eq. (5) total cross sections. This task is carried out in the next Section. 



2.2 Wilson coefficients and inverse Mellin transform 



Resuming all diagrams contributing to leading order in g 2 to the $ - h in- 
teraction, Peskin made explicit the heavy-quarkonium Wilson coefficients 0|. 
They are given by' 

d^-^e^Wr* - 1 ri\<P) 
iV c a (H a + e) 

= 7nlwv\£ ^ii*> (7) 

iV c J (27t) o a [H a + e) 

where \4>) and k are respectively the QQ internal wavefunction and momenta, 
while H s (H a ) is the internal Hamiltonian describing the heavy-quarkonium 
state in a color-singlet (color-adjoint) state, 

k 2 

H s , a = + V s>a (r), 

m Q 

nig being the heavy-quark mass and V s>a the heavy-quark potential. The 
heavy-quarkonium wave function ip(r) in coordinate space and the binding 
energy e appearing in Eq. (7) are determined solving the Schrodinger equa- 
tion 

H S \<P) = -e\<P) (8) 

in the color singlet potential. 



T Note that the coefficients (7) are a factor (e/eo) 2k ~ 1 smaller than in Ref. 
This difference is because the energy A is normalized to the binding energy e in the 
amplitude (3) and not to the Rydberg energy eo as in 
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2.2.1 Coulomb potential 

The leading-twist amplitude (1) was determined assuming the QQ binding 
potential is well approximated by the one-gluon exchange Coulomb potential 



9 2 K 
8 iir 

v a = o (n- 1 ) 



Svrr V c > ' (9) 



in SU(iV c ) gauge theory. To leading order in O \ N C H a is given by the 
free-particle Hamiltonian and the Wilson coefficients (7) read 

16tt f d*k 1 r (f/^g^ 
d ^ = T71 7^3 o I — M ( fc ) 77 „ , / l2 fe-i » 



J (27r) d 3 'fl [{ka o y + e/e ] 

where we have introduced the Rydberg energy e for the system 



/Viv c \ 2 f 

\ 16 7r y m Q aQ 

Solving the Schrodinger equation (8) gives the well-known IS and 2S Coulomb 
wave functions with the corresponding binding energies, 

al ,2 i) {ls \r) = —= exp (-—) ; e ia = e 

V 71 V a / 

which eventually allows for the computation of the Wilson coefficients H 



4 15) =/ 1 -^^ 5/2 (l-^) 3/2 
Jo x 3N; 



(10) 



o 



From Eqs. (6) and (10), the expression for the inverse Mellin transform d(x) 
is straightforward and one gets directly [l(j 

. . 16V 2 (w/e ls -l) 3/2 fl/ 

= ^7T-o \- e J 0(u-e ls ), (11) 



for IS states and 
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ct $( 2S) (u;J - lb x — — a — — ^ 0(a; - e 2S ) (12) 
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for 2S states. Note that these expressions were also obtained by Kim, Lee, 
Oh and Song from the QCD factorization property combined with the Bethe- 
Salpeter amplitude for the heavy-quark bound state, which allowed them to 
include relativistic and next-to-leading order corrections |15j |. 



2.2.2 Screened Coulomb potential 

As stressed in the Introduction, the above formulas may serve as an important 
input to estimate the heavy- quarkonium dissociation process $ + g — > Q + Q 
(or the detailed balance process) in a hot gluon or pion gas formed in high 
energy heavy-ion collisions. We would like here to go one step further and to 
discuss possible medium modifications to these total cross sections. Medium 
effects will be modeled at the level of the heavy-quarkonium potential by 
considering a screened Coulomb potential (Yukawa type) characterized by a 
dimensionless screening parameter //, 



V s = ~T, exp -/ir/a , ( , 

Qur (16) 

V a = 0. 

Solving the Schrodinger equation (8) using the potential (13), the wave func- 
tions and binding energies for IS and 2S states are determined and the corre- 
sponding Wilson coefficients (7) are computed numerically subsequently. For 
the illustration, we plot in Figure 1 the typical size (mean and root mean 
square radii, top) as well as the binding energy (bottom) for the IS (left) and 
2S (right) $ states. Finally, the inverse Mellin transform 

_ ] rc+ioo 

d(x) = — ; — / dz x~ z d(z) 
2 in Jc-ioo 

c being a real constant, is performed thus giving access to the medium-modified 
total cross sections. 



3 Results 



3.1 Finite screening 



Before discussing the results, both the Bohr radius a and the Rydberg energy 
eo in the charmonium and bottomonium channel need to be fixed. Assuming 
both the IS and 2S states to be Coulombic, the heavy quark mass itiq and the 
Rydberg energy eo can be determined from the IS and 2S heavy-quarkonium 



masses. One then obtains 14 
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Fig. 1. Top: Mean radius (solid) and root mean square radius (dashed) of the IS 
(left) and 2S (right) heavy-quarkonium states as a function of jx. Bottom: IS (left) 
and 2S (right) binding energy as a function of ji. 

e 0c = 0.78 GeV , a^ 1 = 1.23 GeV, 
e ob = 0.75 GeV , a^ 1 = 1.96 GeV. 

Using the above (not too hard) scales, the IS (top) and 2S (bottom) heavy- 
quarkonium gluon dissociation cross sections are computed in Figure 2 as a 
function of the gluon energy uj for various values of the screening parameter 
ji. The dominant effect of the screened heavy quark potential is the decrease 
of the IS (respectively, 2S) heavy-quarkonium binding energy from eo (respec- 
tively, Co/4) to e which leads to a lower threshold for the inelastic process. The 
medium modifications of the $ - gluon total cross sections are nevertheless 
not only due to the smaller binding energy, yet the characteristic shapes of the 
cross sections are reminiscent to what is already known for pure Coulombic 
states, fj, = (Figure 2, solid). We checked for instance that the Wilson coeffi- 
cients get somehow modified at finite screening and consequently the partonic 
cross sections do not simply scale as uo/t in Eq. (11). This is a strong indi- 
cation that cross sections cannot be deduced with a simple rescaling of the 
binding energy from 6q to e to mimic medium effects in the heavy-quarkonium 
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dissociation process. Finally, the significant increase of the IS partonic cross 
sections at large screening is particularly noticeable as the dipole size gets 
larger. However, as discussed later, reliable calculations require the space-time 
scales to remain small which prevent one from taking arbitrarily large screen- 
ing parameter values, at least when considering such "light" heavy quarks * . 
Moreover, since the heavy-quark potential in the original QCD analysis needs 
to be Coulomb-like, the screening parameter [i in the model Eq. (13) should 
remain small as compared to one. 
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Fig. 2. IS (top) and 2S (bottom) charmonium (left) and bottomonium (right) gluon 
total cross section as a function of the gluon energy to for various values of the 
screening parameter /i. 

Let us now discuss the heavy-quarkonium hadron cross section. Since heavy- 
quarkonia plunged into the hot medium are most likely to interact with pions, 

* According to ^3], the assumption of heavy-quark Coulombic bound states should 
be appropriate for more than 25 GeV heavy quark mass. 
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we shall only consider the $ - it channel and choose the GRV LO parameter- 
ization for the gluon distribution § in the pion 1^|. The J/ip — 7r and T - tt 
cross sections are computed in Figure 3 as a function of the pion energy A. 
Again, the threshold for the process, located at A = e, gets shifted to lower val- 
ues leading to a strong modification of the heavy-quarkonium pion interaction 
in this region. At high energy, small x = 0(e/X) gluons dissociate heavy- 
quarkonia, thereby increasing the $ - n cross section by a factor (eo/e) 5 where 
5 ~ 0.3 governs the rise of the gluon distribution at small x, xG(x) oc x~ s [l7| . 
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Fig. 3. J/tp — 7r (/e/t) and T — 7r {right) total cross section as a function of the pion 
energy A for various values of the screening parameter [i. 



3.2 Finite temperature 



The $ interaction with gluons and pions has been computed so far using a 
heavy-quark screened Coulomb potential characterized by one parameter fi. 
Interpreting \x as the screening mass in a gluon plasma, the model for the 
finite temperature QQ potential now looks like 

g 2 (r,T)N c 

V s = ex P {-rn D (T) r) (14) 

At short distance and/or low temperature, we shall consider a frozen coupling 
constant 

g 2 {r,T) = g 2 for rT« A (15) 

§ These should be evaluated at a factorization scale e. We take in the following a 
frozen scale eo- 
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and recover the Coulomb potential behavior (9), while the QCD coupling 
starts to run with T at large distance and/or high temperature. At two loops, 
we have 



g\r,T) = ~g\T) 
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In 



2nT\ 51 
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88tt 2 



In 



2 In 



2ttT 



A 



MS / J 



for rT > A 
(16) 

with T c /A]y[g = 1.14 jl8(. The Debye mass m D is related to the temperature 
through the leading-order perturbative result, 



m D (T) = ~g(T) T. 

The A dimensionless parameter introduced in Ref. ^| separates somewhat 
arbitrarily the short from the long distance physics at finite temperature. Fit- 
ting pure gauge SU(3) heavy quark potential, they obtained the empirical 
value A = 0.48 fm x T c . Following we shall take the 2-loop running cou- 
pling (16) rescaled by 2.095 and interpolate smoothly between the short and 
long distance regime ' . 

The partonic and hadronic J / ip and T cross sections are computed in Figure 4 
for several temperatures in units of the critical temperature for deconfinement, 
T c = 270 MeV in SU(3) pure gauge theory |2j. The temperatures selected for 
the bottomonium system are chosen to be slightly higher than those for the 
charmonium system since the larger bottom quark mass (hence, smaller size) 
probes more efficiently hotter QCD media [19j. 

The effects of the running coupling in Eq. (14) being quite small, rather similar 
features at finite temperature and at finite screening are observed. In partic- 
ular, the charmonium binding energy (hence the inelastic threshold) drops by 
a factor of two already at T/T c = 0.5 and thus affects dramatically the J/ip 
interaction in the vicinity of the threshold. At higher temperature, the J/ip - 
gluon cross section is significantly enhanced at small gluon energy due to the 
larger charmonium size. The J/tp — ir cross section is also somewhat modified 
with a magnitude increasing noticeably with the temperature. Moving to the 
bottom sector (Figure 4, right), the T cross sections exhibit the same general 
characteristics yet the medium effects at a given temperature prove much less 
pronounced from the smaller bottomonium size. 

At high temperature, heavy-quarkonium interaction can not be described by 
short-distance techniques (see Fig. 1) and our predictions are not valid any 
longer. On top of that, the process described here is the heavy-quarkonium 
dissociation by hard gluons as opposed to the soft gluons which only af- 
fect its properties. Therefore, our calculations should be valid as long as 



™ Similar results are obtained using the one loop running coupling with an appro- 
priate rescaling. 
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Fig. 4. J/tp (left) and T (right) total cross sections with gluons (top) and pions 
(bottom) at various temperatures. 

the Debye mass is kept smaller than the heavy-quarkonium Rydberg energy, 
m D(T) < eo- This condition is fulfilled provided the bath temperature is 
smaller than 350 MeV. Above that scale, the screened exchanges are able to 
dissociate the bound states, the factorization between the heavy-quarkonium 
physics and the external gluon field is broken and the above QCD picture loses 
its significance. 



3.3 J /if) spectral function width 



The former results indicate that Debye screening effects may play an important 
role in the heavy-quarkonium dissociation by incoming gluons or pions. In 
order to illustrate how medium modifications could affect the $ suppression 
in heavy-ion collisions, we compute in this section the IS charmonium thermal 
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width r j/Vi (or equivalently its lifetime, r Jfii = V * ) in a hot gluon bath. 
Assuming the J /if) suppression is only due to the gluon dissociation process, 
the width can be written 

1 f°° 

where n g {uj,T) = 2(N% — 1) j (exp {oj/T) — 1) is the gluon density in a gluon 
gas in thermal equilibrium. 



^ 0.3 
> 

0) 




T/T c 

Fig. 5. J/i/> thermal width as a function of the temperature with {dotted) and without 
{solid) modifications of the heavy-quark potential. The lattice data point obtained 
in Ref. |2Q] is also shown for comparison. 

The thermal width is computed in Figure 5 as a function of the tempera- 
ture T assuming the vacuum {solid) and the in-medium {dashed) J/ip - gluon 
cross section. At small temperature, T<6, most gluons are not sufficiently 
energetic to dissociate J/ip states and the width remain small as the phase 
space selected by the J/tp gluon threshold is restricted. When the medium 
gets warmer, more and more gluons are able to interact inelastically with the 
J/ip, hence the thermal width increases. Interestingly enough, the in-medium 
J I ip thermal width proves larger by a factor of two or more up to T = T c due 
to the lower threshold in the medium modified cross sections. At even higher 
temperature, the medium modified result becomes smaller to that in the vac- 
uum since dissociating gluons (with uj of order e) grow scarce. Also plotted in 
Figure 5 is the J/ip width computed recently on the lattice at finite temper- 
ature in the quenched approximation j^. Although a significant discrepancy 
remains between our calculations and the lattice data point, it is interesting 
to note that adding medium effects tends to reduce the disagreement, whose 
origin is not clarified. 
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4 Concluding discussion 



Before summarizing our main results, we would like to discuss the limitations 
of our approach. The starting point of the calculation is the forward scattering 
amplitude Ai<s>h originally derived for Coulomb bound states. To go beyond 
this one-gluon exchange picture would require to include light quark loops 
in the diagrammatics, to which the soft gluon source may couple, that we 
have not attempted. However, as conjectured in it is appealing to guess 
that the generic dipole coupling appearing to leading order in g 2 in the heavy- 
quarkonium Wilson coefficients (7) survives perturbative and non-perturbative 
modifications of the QQ binding potential. Therefore we believe that taking 
the literal expression for the Coulomb states Wilson coefficients and compute 
them in a screened Coulomb potential appears sensible, at least as long as 
the screening remains reasonable, tudCLq <C 1. This is certainly the case when 
the temperature is kept small as compared to the heavy quark mass. In that 
sense, the smallness of the charm and bottom quark mass as compared to 
the non-perturbative scale of QCD indeed remains a problematic issue. As 
we have seen, typical space time scale becomes increasingly larger with the 
temperature, thus strongly limiting our confidence in the high temperature 
regime. Finally, one should keep a clear factorization between the gluon source 
and the heavy- quarkonium swimming in the gluon bath. We have seen that 
such a separation should be achieved as long as the Debye mass is small 
as compared to the bound state Rydberg energy, that is for temperatures 
T < 350 MeV. 

We presented a numerical calculation of the heavy- quarkonium cross section 
with gluons and pions, taking into account the possible medium-modifications 
of the heavy-quark potential at finite temperature. Such a work can therefore 
be useful to estimate heavy-quarkonium production in high energy heavy- 
ion collisions. In particular, we feel it would be interesting to explore the 
phenomenological consequences of such corrections comparing them to present 
calculations based on the vacuum heavy-quarkonium interaction. Finally, this 
very framework could be applied to study the $ interaction using a variety of 
realistic heavy quark (confining) potentials currently used in charmonium and 
bottomonium spectroscopy [l2| to describe more accurately, although further 
away from the perturbative requirement, heavy-quarkonium interaction with 
gluons and hadrons. 
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